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Abstract 

A non-linear sigma model mimicking the Skyrme model on 53 is 
proposed and a family of classical solutions to the equations are con- 
structed numerically. The solutions terminate into catastrophy-like 
spikes at critical values of the Skyrme coupling constant and, when 
this constant is zero, they coincide with the series of Harmonic maps 
on 53 constructed some years ago by P. Bizon. 



1 Introduction 

A few years ago P. Bizon [|IJ constructed a remarkable sequence of harmonic 
maps from S3 into S3 . Only a few examples are known of non-linear equations 
which admit such a sequence of regular solutions. Some of them occur in the 
framework of gauged non-linear sigma models. The prototype of non-linear 
sigma model in theoretical physics in the Skyrme model . It is based on a 
mapping of the physical space R 3 into the manifold of the Lie group SU(2), 
wich is homeomorphic to S3. In order to guarantee the existence of topolog- 
ical soliton the physical space is compactified. The original Skyrme model 
possess an SU(2)®SU(2) global symmetry which can partly be gauged. Sev- 
eral examples of sequence of regular-finite energy classical solutions appear 
in different gauged vesions of the Skyrme model f||, [§]. 
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For all known examples, an interesting phenomenon, related to catastro- 
phe theory, occurs to the solutions when the minimal (i.e. quadratic in the 
derivatives) sigma model is supplemented by a quartic term in the derivatives 
[@> H) 0) a so called Skyrme term. The higher derivative interaction term is 
proportional to a free constant, say e, which constitutes a coupling constant 
of the model. To be more concrete let us denote (in the models |5|, |], ||) 
by S n (x, e — 0), n G IN the sequence of solutions available in absence of the 
Skyrme term and E n (0) the corresponding classical energy. Solving the clas- 
sical equations for e > 0, it appears that a branch of solutions S n (x, e) exists 
for < e < e cr (n) and that no solution is available for e > e cr (n). However, 
there exists a second branch of solution, say S n (e), for e G]0,e cr (n)]. The 
numerical results indictate that the following properties hold for n G N 

E n {e) < E n {e) (1) 

lim E n (e) = E n (e) (2) 

e— >e C r(n) 

lim£ n (e) = £ n _i(e) (3) 

e— >0 

In fact the statement (2) is even stronger because we have in fact 

lim S n (x, e) = lim S n (x, e) = S n (x, e cr (n)) (4) 

in a uniform sense for x G R 3 . The counterpart of Eq.(4) related to the 
limit (3) does not hold because the solutions S n (x, e) and S' n _i(x, e) do not 
in general obey the same boundary conditions. 

The plot of the energies E(n, e) and E(n, e) as a function of e shows 
that the two curves terminate at e = e cr (n) into a cusp, characteristic of 
catastrophe theory. 

This phenomenon is absent in the normal (i.e. ungauged) Skyrme model 
but it occurs in several gauged versions of this model, namely in the "hidden 
gauge" Skyrme model [0. In passing let us point out that the trivial topo- 
logical sector sector of the "hidden gauge" Skyrme model is equivalent to the 
(bosonic part of the) electroweak standard model considered in the limit of 
an infinitely heavy Higgs particule (i.e. for Mh ~ A — ► oo) |5|, ||. 

Owing the above properties associated with sequences of classical solu- 
tions available in the gauge-Skyrme model, it is natural to study how the 
harmonic maps of Ref I] respond to the addition of a Skyrme interaction 
term. This is the aim of this paper. 
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2 The model 



The functional energy studied in reads 

where / represents a mapping from S3 (with local coordinates noted x l ) into 
S 3 (with local coordinates note f A ). We use the same notations as in [p]], 
namely x l = (if), 9, (p) are the standard coordinates on S 3 with the metric 

ds 2 = dip 2 + sin 2 ip(d8 2 + sin 2 6d(p 2 ) (6) 

and similarly f A = (\I>, 6, $) on the S 3 target space. 

The functional Eo(f) is a non-linear sigma model; the classical solutions, 
i.e. the extremum of the functional, are called harmonic maps in mathemat- 
ics. In this paper, we study the extremum of the extended energy functional 

E(f) = E (f) + E Sk (f) (7) 

with the Skyrme interaction term 

E M ) = I hMf)hcM) ^ a £- d £ { , ig ^ + , 2gW (8) 

where fix, \i<i are constants. In order to keep the model (H) as close as possible 
to the standard Skyrme model defined on IR 3 we choose jii = — /i 2 = k 2 . 

Finally, we also assume that the function f A (x % ) is equivariant i.e. in- 
variant under the SO (3) rotations corresponding to the angles 9 and (p. This 
leads us to the ansatz 

*(V0 = /(V0 , 6 = , $ = (9) 
in terms of which the functional E(f) reads 

£(/ =4vr / {sin ip{—) +2— ^— + /? (— -g— (— ) + 4 )}# (10) 
jo sin -0 sm ip dip 2 sin -0 

The classical equations are obtained by varying E(f) with respect to /. We 
are interested by solutions such that the energy density (i.e. the term in 
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bracket in the integral above) is regular. This imposes in particular /(0) = 
/(tt) = 0. 

In the purpose of comparison with the normal Skyrme model on R 3 we 
just write the expression of the classical energy in this model: 

E(f) = 4, f *{(f? + 2^ + .^(ff + ^-)) i r (11) 
jo ar r z r z ar Zr^ 

3 The numerical solutions 

It is convenient to label the extrema of the functional (|T0D by f n (ip, k 2 ), where 
the integer n counts the number of times the solution crosses the line / = f . 
It is also useful to notice the discrete invariances of the equation : 

/(V0-/(tt-V0 , /(V)-^7r-/(7r-V) • (12) 

They transform any solution into another one with the same energy. 

In a sequence of extrema f n (ip,0) (also called critical maps) to the 
functional (||) was constructed. It is such that fo(ip, 0) = , fi(ip, 0) = if). To 
our knowledge, the solutions corresponding to n > 1 cannot be expressed in 
terms of known functions. The profiles of the functions f2{fp, 0) and f? J {4 ! i 0) 
are represented on Fig.l and Fig. 2 respectively. For the first few values of n, 
the energies are given || in the table below : 



n 


E(n)/6rr 2 








1 


1 


2 


1.2319 


3 


1.3036 


4 


1.3240 



We now discuss the evolution of these solutions for k 2 ^ 0. Obviously the 
vacuum solution f = (or, by fll2]) , / = it) exits for all values of k 2 and its 
energy is zero. The solution corresponding to n = 1 reads 

A(^ 2 ) = V> or / 1 ^, K 2 ) = 7r-V (13) 



4 



and, like the vacuum, also exists irrespectively of the constant k 2 . It has an 
energy 

£ 1 (« 2 ) = 67r 2 (l + ^) (14) 

Solving numerically the equations for finite values of k 2 we find that 
the solution f n (i[),0) gets continuously deformed and generates a branch of 
solutions f n (ip, k), at least for small values of k. For n = 2 the solution exists 
on the interval 

< k 2 < k 2 ct {2) ; k 2 ct {2) ~ 0.261 (15) 

E 2 (K 2 cr (2))/Qir 2 ~ 1.670 (16) 

The profile of f 2 (ip, K l r {^)) * s a ^ so shown on Fig.l and the energy of the 
solution is represented (in function of k) by the upper part of the n = 2 
curve in Fig. 3. We find no solution for k, 2 > K 2 r but a second branch of 
solutions, say f 2 (ip, ^ 2 ); exists on the interval (|15|). For fixed k, the energy of 
the solution f 2 is lower than the energy of f 2 , i.e. (using an obvious notation) 

E 2 {k) < E 2 (k) (17) 

as represented by the lower and upper branches of the n = 2 curve in Fig.3. 
At k 2 = «cr(2) the solutions f 2 and f 2 coincide and the energy plot terminates 
in a cusp shape, suggesting some catastrophic behaviour. 

We would like to stress that the evolution of the solution profile from 
the upper branch to the lowest one is completely smooth, the profiles of the 
solution f 2 (ip, 0.2), f 2 (ip, 0.002) are displayed of Fig.l. The numerical analysis 
strongly suggests (as illustrated by Fig.l) that 

lim f 2 (ip, k 2 ) = ix for < ip < tx (18) 

This statement is also supported by the evaluation of the energy 

lim E 2 (k 2 ) = E (0) = 0. (19) 

K 2 — >0 

as indicated on Fig.3 

Similarly, the branch developping from the n = 3 solution exists on the 
interval 

< n 2 < K 2 cr (3) ; K 2 cr (3) ~ 0.017 (20) 
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and 

£3(4,(3))/6tt 2 ~ 1.395 (21) 

Again the second branch of solutions f^ip^n 2 ) occurs and a few of these 
solutions are superposed on Fig. 2 (resp. fsfy, 0), fs(ip, « 2 r (3)), fzi^i 0.015), 
fs(il>, 0.0005)). The numerical results suggest that 

lim f 3 (ip, k 2 ) = 7r - V for < ip < vr (22) 

to be compared with (pT3|). We also find lim K 2^ E 3 (k 2 ) = Ei(0) = 6n 2 . 

In view of these results, we conjecture that two branches of solutions say, 
/ n (/€ 2 ) and f n (n 2 ), exist for < k 2 < ft 2 r (n) and that the result suggested by 
our numerical analysis holds for n > 3, that is to say 

lim vUk 2 ) = / n _ 2 (0) for < ^ < vr (23) 



K->0 



lim E n (n 2 ) = E n _ 2 (0) (24) 



4 Stability 

The stability of the solutions constructed in the previous section can be 
studied by perturbating the classical solution, say f(ip)- In this paper, we 
limit ourselve to fluctuations depending only on ip, i.e. 

m-iw+vw (25) 

Expanding the functional (10) in powers of r/ leads to a quadratic form which 
can be diagonalized. The relevant eigenvalue equation reads 

2 sin2 ./\ d 2 



2j/ sin2 ./\ d 



sin 2 ip dip 2 

-'•<>"• ''•'"'/ . , 

sin %p dip 

COS 2/ , 2,Fr,2 COs2 f 



+ (2-^L + K 2 (f 

\ sin 7/j 



sin" 1 7^ yJ ' sin 2 -0 



+ ^-^-(sin 2 2/ + 2 sin 2 /cos 2/)) V = cu 2 V (26) 

/* si n in ' 



2 sin ip 
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and determines the normal modes r](ip),u) 2 about the classical solution f(i/j). 

For the vacuum, f(if>) = 0, this equation becomes independent of k 2 and 
can be solved explitely; the normal modes read 

rj k (x) = sin^ Cf(cos?/>) , w£ = 3 + k(k + 4), (27) 

with k = 0, 1, 2, . . . and where Cf(t) denote the Gegenbauer polynomials 

C 2 = l , C 2 (t)=t, (28) 

(n + 2)C 2 +2 (t) = 2(n + 3)tC 2 +1 (t) - (n + 4)C 2 (t) (29) 

The lowest eigenvalues associated with the solutions f(ip) = f n (ip,0) are 
given in [0 . Apart from the case n — 1 whose eigenvector have the form (p7|) 
with eigenvalues 

cj 2 = -1 + k{k + 4) (30) 

the spectrum has to be computed numerically. In particular the solution 
f n (ip,0) possesses n independent modes of instability of the form (p5|). 

We studied Eq. fl26|) numerically for k 2 > and for the two modes 
of lowest eigenvalues associated with the solutions corresponding to n — 1 
and n = 2. Our results are reported on Fig. 4. The solid lines a, b (resp. 
c,d) represent the modes (P7j) , (resp. (|3"0"D) for = 0, 1; we see that the 
negative mode about fi(ip,K 2 ) becomes positive for k > 0.8 (this is not 
sufficient, however, to claim that fi(ipiK 2 ) is a local minimum of (7), general 
fluctuations have to be considered as well). We also notice that the solution 
/2(Vs ft 2 ) has no directions of instability in the channel (^5|). 

The two lowest modes relative to the solution ^(V*; ^ 2 ) (resp. /^(Vs ^ 2 )) 
are represented by the dashed (resp. dotted) lines on Fig. 4. The numerical 
results strongly support the fact that the lowest modes about the solution 
f'2(ip,K 2 ) and the solutions /^(V^ 2 ) vanish in the limit k 2 — ► ^(2) This 
confirm on an explicit example a general expectation (based on catastrophe 
theory; see. e.g. 0) that the solution on the higher branch (i.e. f%) has one 
more directions of instability than the solution on the lower branch (i.e. ^2)- 

5 Conclusions 

The non-linear sigma model based on the mapping of S 3 into S 3 admit a 
rich set of classical solutions [|TJ . By supplementing it by a Skyrme-like term, 
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we have shown that the solutions can be interpolated by the extrema of the 
family of models indexed by the coupling constant k 2 . 

The energy of the classical solution then leads to a serrated pattern sum- 
marized by Fig. 3. Apart from the two lowest branches, which respectively 
play the role of the vacuum (E = 0) and of the Skyrmion solution of stan- 
dard Skyrme model, the other solutions cease to exist at critical values of the 
parameter k 2 , where two branches of solutions terminate into a catastrophic 
cusp. 

The analysis of the normal modes about the different solutions indicates 
that their stability strongly depends on the Skyrme term and suggests that 
the spectrum of the quadratic fluctuation operator varies continuously with 
the parameter k 2 . 

Acknowledgements We gratefully acknowledge discussions with P. Bi- 
zon. 

Figure Captions 

Figure 1 The solution faty, k 2 )/tt is plotted in function of ip/ir for four different 
values of the coupling constant k. 

Figure 2 The solution fz(ip, k 2 )/7t is plotted in function of tp/n for four different 
values of the coupling constant k. 

Figure 3 The energy of the solutions corresponding to n = 1, 2, 3 is plotted in 
function of the coupling constant k. 

Figure 4 Some eigenvalues of Eq.(p6|) are plotted in function of k. The lines a, b 
refer to the vacuum (p7|) and the lines c, d to the n — 1 solution. The 
dashed (resp. dotted) lines refer to the modes of the solution / 2 (resp. 

h). 
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